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Abstract. Unified description on the long-time tail of velocity autocorrelation function (VACF) and the 
algebraic decays of the equal-time spatial correlation functions for nearly elastic and uniformly sheared 
granular liquids is developed based on the generalized fluctuating hydrodynamics. We predict that the 
cross-over of the long-time tail of VACF from t~ 3 ^ 2 to t~ 5 ^ 2 with the time t regardless of the density. We 
also demonstrate the existence of algebraic tails of the equal-time spatial density correlation function and 
the equal-time spatial velocity correlation function which respectively satisfy r~ 11/ ' 3 and r~ 5//3 for large 
distance r. 



PACS. 61.20.Lc Time-dependent properties; relaxation 
- 45.70.-n Granular systems 



05.20.Jj Statistical mechanics of classical fluid 



1 Introduction 

Many theoretical aspects of granular assemblies have been 
described by gas kinetic theory such as Boltzmann-Enskog 
equation [TJ[2] in which effect of correlations only appears 
through the radial distribution function at contact. Based 

, on Chapman-Enskog method or Grad expansion, it is pos- 
sible to determine the transport coefficients and to derive 
hydrodynamic equations 3 ( 4 , S1E] - I n spite of the success 
of semi-quantitative description of granular assemblies by 

' the kinetic theory [7] , recently we have recognized impor- 
tant roles of long-time tails and long-range correlation 
functions in granular liquids [51[TJ|[TOlP1[T^[n|PllP31[Trjl 
[TTlfTS] which cannot be described by Boltzmann-Enskog 
theory. 

In general, correlation effects are complicated, which 
, strongly depend on the boundary condition. However, we 
can extract simple statistical properties of uniformly sheared 
granular liquids which are much simpler than nonuniformly 
sheared systems and independent of the boundary condi- 
tion. Indeed, recent analysis on uniformly sheared granular 
liquids with small inelasticity has revealed the existence 
of long-range spatial correlation function, the suppress 
of the long-time tail of velocity autocorrelation function, 
the existence of generalized Green-Kubo formula, and the 
integral fluctuation theorem P^IHTIH^] . This uniformly 
sheared granular system can be achieved, at least, even for 
boundary-driven dense flow of frictionless granular parti- 
cles [20j. Moreover, we need to remove the boundary ef- 
fects to extract the spatial correlation function of granu- 
lar liquids which should be translational invariant. How- 



ever, our previous analysis on the long-time tails assume 
that granular gas is dilute, while experiments analyze the 
dense flows. So far, the connection between the long-time 
tails for velocity autocorrelation function and the equal- 
time spatial correlation functions is not clear. Hence, we 
demonstrate both the long-time tails and the long-range 
correlation functions for uniformly sheared granular liq- 
uids with small inelasicity can be derived by an unified 
method based on the generalized fluctuating hydrodynam- 
ics gH|31231|23|. 



In this paper, thus, we theoretically investigate the 
velocity autocorrelation function (VACF) and the equal- 
time correlation functions based on the generalized fluctu- 
ating hydrodynamics with a linear non-local but instante- 
neous constitutive relation and fluctuation-dissipation re- 
lation. As stressed in Refs. 1 .''>. 1 7 2"> . the formulation can 
be used for any sheared isothermal liquids besides granular 
liquids. However, to clarify the argument, we only focus on 
uniformly sheared granular liquids with small inelasticity. 
In Sec. [H we will summarize the outline of the general- 
ized fluctuating hydrodynamics. In Sec. [3] we will analyze 
the linearized generalized fluctuating hydrodynamics. In 
SecSl we will calculate the correlation functions and their 
asymptotic forms in the long-time limit and the large dis- 
tance limit. In Sec. [5] we will discuss and conclude our 
results. In Appendices, we present some details of our cal- 
culations. 
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2 Generalized fluctuating hydrodynamics 

The systems we consider are three-dimensional sheared 
granular flows consisting of TV identical smooth and in- 
elastic hard spherical particles in the volume V, where 
the mass and the diameter of each grain are respectively 
given by m and a. We should note that granular systems 
often cause shear-bands near the boundary. Since we are 
not interested in complicated behaviors under such shear 
bands, we focus on the uniformly sheared granular liq- 
uids characterized by a component of the velocity field 
c a (f) = iy$a,x with the shear rate 7. It is known that 
such uniform flow can be realized if the system size is not 
extremely large and granular assemblies are nearly elastic 
under Lees-Edwards boundary condition. We also assume 
that the particles collide instantaneously with each other 
by a restitution constant e which is less than unity. 

Let us consider the velocity autocorrelation function 
in the three-dimensional sheared granular liquids, defined 
by 



Since we assume uniformly sheared case, the tempera- 
ture T only appears through the time evolution equation 
of the velocity fields u(r,t). Thus, the equations for the 
generalized hydrodynamic equations are given by [31,22, 

d t n + V • (nu) = 0, (7) 



1 



1 



d t u a + u a Vpu l3 + — VaM + — V,3 (SSb + ZSb) =0, (8) 

where n is the number density. Here, we have introduced 
the generalized chemical potential or the effective pressure 



/' 



T 



In n 



dr'Cir - r', e, j)5n(r' , t) + 



, (9) 



where C(r — r',e,7) is the direct correlation function, 
which satisfies nC(k,e, A f) = 1 — S(k,e, A f)^ 1 with the 
structure factor S(k, e,j). E^(r,t) is the viscous stress 
tensor given by 



1 N 

C S ,aa(t) = —^2(5Vi ta (t)5Vi, a (0)), 



(1) 



where Svi(t) = Vi(t) — c(ri(t)) is the peculiar velocity of 
the particle i. In addition, we also consider the equal-time 
spatial correlation functions under the steady state 

C nn {r) = (5n(r + r',0)Sn(r',0)), (2) 
C pp (r) = {Sp(r + r',0)-6p(r',0)) (3) 

for the fluctuations of the number density Sn(r, t) = n(r, t)- 
uq and the momentum density Sp(r, t) = m n(r, t)5u(r, t) 
with the fluctuation of the velocity field u(r,t) given as 
Su(r,t) — u(r,t) — c(r). 

Introducing the Fourier transform of any function f(q) = 
J drf{r)e lq r , and following Ref. |2B], we can rewrite Cs, aa (t) 
as 



SSp(r, t) = ~J dr' [ V (r - r', e)e a0 (r', t) 

-{2 V (r - r', e)/3 - C(r - r', e)}e 77 (r', t^ftO) 

where e a p(r,t) = {V Q u (3 (r,t) + Vpu a {r, t)}/2. S^ (r,t) 
is the random part of the stress tensor satisfying (E^p) = 
0, and the fluctuation-dissipation relation (FDR) 



EUr,t)Z* 5 (r',f)> = 2T6(t - t'){rj(r - r',e)A a 



+C(r - r',e)S al3 6 7 s} 



0fS 

(11) 



with A a ps<y = Sa-ydps + 5 a s5p 7 — 25 a p5 1 s/'i- Here, we have 
used Einstein's sum rule on the Greek subscript. The gen- 
eralized shear viscosity rj(r, e) and the generalized bulk 
viscosity C( r : e ) are represented by v*(k,e), v^ik.e) and 
Enskog's mean free time 



t E 



1 



/ TY17T \ 

4irnoa 2 go((T,e) V T ) 



1/2 



(12) 



C s , aa (t) = J -^-^{5u a {q, t)Su a (q',0)) 



(2tt) 3 (2tt) s 
(P(q,t)P(q',0)), 



with the radial distribution function go(a, e) at contact as 
f*(k, e) = (mn cr 2 i i j 1 )~ :L (C(fc, e)+4^(fc, e)/3), and v\ (k, e) 
(4) (mno<r 2 t] s 1 )~ 1 ri(k, e), where no, T](k, e) and C(fc, e) are re- 



where P(q,t) is the Fourier transform of the the micro- 
scopic concentration of the tagged particle i defined by 
P(r,t) = S(ri(t) - r) (see Appendix E}. 

Similarly, the spatial correlations in Eqs. © and 
can be rewritten as 



spectively the average number density, and Fourier trans- 
forms of 7y(r, e) and £(r, e). 

It is known that v*{k,e) and v%(k, e) are respectively 
given by 



1) = 2(1 - j (k) + 2j 2 (fc))/(3fc 2 ), (13) 



dq dq' 

(2tt) 3 (2tt) 3 



(8n(q,0)Sn(q',0)) 



and 



^(fc,l) = 2(l-j (fc)-j 2 (fc))/(3fc 2 



(14) 



-iqr 



(5) 



c pp( r ) - ( m n o) 2 



dq dq' 

(2tt) 3 (2tt) 3 



(Su(q,0) ■ Su(q',0)) 



^ — iq-r — i(q+q')r' 



(6) 



where we have used the Fourier transform of 5n(r, t) and 
Sp(r,t), and 8p(r,t) ~ m noSu(r,i). 



for elastic hard spherical particles, where ji{k) with / = 
or 2 is the Z-th. order spherical Bessel function [2T1[2^1 
l25 P ^ | l27 p g ] . Although we do not know how v{(k, e) and 
1^2 (k, e) depend on e, the explicit e-dependences of vl(k, e) 
and e) are not important in this paper. Therefore, 
we are keeping discussion without their explicit forms. 

It should be noted that the equilibrium or the un- 
sheared structure factor So(k, e) = S(k, e, 7 = 0) is given 
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by an approximate expression of the pair-correlation func- 
tion for unsheared granular liquids |29j . which covers the 
equilibrium pair-correlation in the elastic limit. In Eq. 
(|10j) . we assume a linear and instantaneous transport law, 
which is valid when the inelasticity of the particles is 
small in the granular liquid. If the inelasticity is not small 
enough, we need to use a non-linear transport laws because 
the normal stress differs from the shear stress |30) . This 
set of equations (|T))- (fTTj) is a reasonable starting point for 
isothermal molecular liquids (e = 1), once we use appro- 
priate generalized transport coefficients and So(k,e = 1), 
and can be used to describe isothermal sheared granular 
fluids in the vicinity of e = 1 [17] . We also note that FDR 
(fTTj) is not satisfied in highly dissipative granular liquids. 
Indeed, we can demonstrate the existence of some correc- 
tions in the expression of the noise correlation for general 
case. It is possible to derive the fluctuating hydrodynam- 
ics from the Liouville equation and its corresponding Mori 
equation with the Markovian approximation and drop- 
ping some correlated noise terms [31] . The use of FDR, 
however, can be justified if the system is nearly elastic 
and weakly sheared situation and the density is far from 
jamming transition point (point J). The validity of the 
generalized fluctuation hydrodynamics has been verified 
through the comparison between the theory and the sim- 
ulation if the density is far from point J [THE]- We ex- 
pect that non-Markovian effects play important roles near 
point J. 



3 Linearized equations around uniform shear 
flow and their solutions 

In this section, we analyze the hydrodynamic fluctuations 
given by Eqs. © and §5§. In the first part, we introduce 
the linearized equations for the fluctuations. In the sec- 
ond part, we explicitly write the solution of the linearized 
equations. 



3.1 Linearized equations 

Let us introduce the non-dimensionalized vector z(r,t), 
whose Fourier transform is given by 

z T (q,t) = (5n(q,t),^ s ( g) t)/(t-V 4 ), 

6u y (q, t)/{t E l cr A ), Su z (q, t) / {t~ E l a 4 )) . (15) 

From Eqs. (TTJ) , (THJ) , and (fl"5")) . we obtain the linearized evo- 
lution equation for z(k,t) = z(q,t) as 



dt-rk,^)z I z H. .I'm 



where the time, the wavenumber and the shear rate have 
been non-dimensionalized by t = tsi, q = fc/c, and 7 = 
7*/tE, respectively. Here, the matrix L is expanded as 



where l_o and l_i are respectively given by 
Lo = — L50 + Llq + Ldo 



with 



no<j 3 ik x no<7 3 iky no<7 3 ik z 
p*ik x 

p*ik y 

p*ik z 



(18) 



(19) 



-LO 



u*k z 
v Z K x 



is^kxky i / ^k x k z 



ky kx ^3 ky ky k z 
Vn k z k x Vn k z ky is*, k 2 



(20) 



-150 — 





Pnk 2 






v*k 2 

-* 1 - 2 



v*k 2 



(21) 



and 



Li 




10 





(22) 



with 1/3 = v\ — 
A = T/(mn a 5 t^ 
ponents 



p*(fc,e) = ASo(*,e) _1 and 
). The random vector R has four com- 



Pi 



Ri = 0, 
R a+1 = (mn aH E ' 



(23) 



where a = 1,2,3 respectively correspond to x,y and z. 
Although 7 dependence of S{k,e, A f) should appear in Lj, 
we simply ignore such terms. The validity of this simpli- 
fication has already been checked from the comparison of 
the results with our simulation (17j . 

We should note that the wave number k a is always 
larger than A — 2-k/L with the system size L, because 
the uniformly sheared flow becomes unstable for sheared 
granular liquids in the limit k — ¥ 0, and we need infinitely 
large work to produce the uniformly shear flow for all sit- 
uations. This singular behavior around k x = can be 
understood because the effect of the shear in the left hand 
side of Eq. (TTB1) disappears if we take k x = 0. Thus, it 
is not appropriate to take the limit k x — ¥ without the 
introduction of the infrared cutoff A, 



3.2 Solution of the linearized equations 

The solution of the linearized equation (|16p is easily ob- 
tained following the parallel procedure in [P71I25] . We in- 
troduce the right eigenvector ip^\k) and the eigenvalue 
A^(fe) satisfying 



L = L +7*L! 



(17) 



_d_ 

dk,. 



ip^(k) = X u \k)ip U) (k). (24) 
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We also introduce the associated biorthogonal vector, i.e. where 

the left eigenvector (p^(k), satisfying ip^'(k) ■ ip^\k) = j 

Skj- E^ j) (k,t} = exp[- dfA^(fc(7*f))], (39) 

Following the procedure described in Appendix iBl we Jo 

obtain the eigenvalues . , 

with k(r) = (k x , k y + rk x , k z ). Hence, the solution of Eq. 

. m . . , k x k y (1) (HH) can be formally represented by 

A — A + + 7 r 2 s l fc Ji l^J 

A (2) = A_+ 7 *^e (2) (fc), (26) *(fc,«)=y;/ rfsr'^t-^H^r^-s)).^, 

A( 3 )^^(fc,e)fc 2 - 7 *^, (27) u ( 4 °) 

Ar where 

A (4) =^(fc,e)fc 2 (28) 



tp U \k,t) ee i/>«(fc).E«(fc,i). (41) 



within the approximation up to 0(7*), where we have 
introduced 



v\(k,e)k 2 + yj(v* 1 {k 1 e)k 2 )' 2 - 4n Q <j 3 p*(k, e)k 2 
X+ — 



4 Correlation functions 



This section is the main part of this paper, in which we 
(29) 

present the explicit forms of correlation functions. This 



_ f*(k, e)k 2 - y/ (v{{k, e)k 2 ) 2 - 4n Q <r 3 p*(k, e)k 2 section consists of three parts. The first part summarizes 

2 the general results of correlation functions. In the second 

(30) part, we evaluate the long-time tail of the velocity au- 
tocorrelation function. In the third part, we discuss the 
£<4)(fc) = A + 4- n o a " K " hQ k p*(h e ) (31) equal-time long-range correlation functions which are es- 

v ' s ~' " sentially the same as that in Ref. [T7] . 



\2 
A + 


no<j 3 k 2 


N 2 


h 2N 2 


X 2 _ 


n a 3 k 2 


W~ 


h 2N 2 _ 



^{k) ee + J^kd k p*(k,e), (32) 



and 



N% = -n a 3 p*(k, e)k z + X z + , (33) 
N 2 _ = -n <r 3 p*(k,e)k 2 + X 2 _. (34) (z a (k, t)zp(k', 0)) 



4.1 General results for correlation functions 

In order to obtain the correlations from Eq. ((U)-©, we 
consider (z a (k, i)zp(k , £)), which is represented by 



1—1 r* r° 

See Appendix IBl for the detailed expressions of the eigen- _ / ^ g ds'^ l \k i—s)$ m \k' —s') 

vectors. , , «/-<» J -00 ° P 7 



l,m—l ' 



Using the eigenvectors ip^\k) and tp^(k), the vector 
z(k,t) is represented by x (F«(fc( 7 *(t- s)), s)F^(k ( 7 *(-*'))> «')>• ( 42 ) 



z(k,i) = J2a U Hk,t)*l>(k), (35) 



FromEqs. ^ and (TSTJ), (FW(fe,l)F( m )(fe',t')> satisfies 



(F (l \k,t)F^ m \k',i')) = {2TT) 3 S 3 (k + k')5(i-i')F^ m \k), 

(43) 

where we have introduced a^(k, t) = <p^'(k) ■ z(k,t). with 

Substituting Eq. (1551 into Eq. ([TH]) and taking the inner ^2 
product with ip^\k), we obtain the time evolution cqua- F^ lx \k) = — 2Ak 2 v\(k, e )^jj 



tion of a^\k, t) as 



9 ™, A ^ ,,,, , . F^(k) = -2Ak 2 ^(k,e 



(dt-rk x ^ + \ (j) (k^j a^(k,t)=F^(k,t), (36) 



2 



2 



F< 12 )(fc) - F< 21 >(M~) = -2Ak 2 vl{k,e)^- 
where we have introduced A + A_ 

FW(M-) - • HflM"). (37) F(33)(fc) = 2Ak ^ k ^ 



As shown in Appendix [C] the solution of Eq. (|3"6")l is 

given by 



F< 44 )(fc) = -(A/(fc) 2 + l)F( 33 )(fe), 
F< 34 )(fc) = -F^ 3 \k) = M(k)F^ 33 \k), (44) 



and 



a^(k,t)= / d Sj B«(fe,t- S )F«(fe( 7 *(f- S )), S ), (38) F« m \k)=0, 



(45) 
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for other pairs of (I, m). Hence, we obtain 
(z a (k,t)z e (k' ,0)) 



J2 f dsf ds'^(k,t-s)^ m \k\- s ') 

i J — oo J — oo 



l.m— 1 ' 



x(27r) 3 ,5 3 (fe(7*(i- 8 ))+k(ff*{-8'))) 
x5(s- s')F^(k(j*(i- s'))) 
= (2n) 3 5 3 Ck(rt) + k')C aP (k,t), (46) 
where we have introduced 

C aP (k,t)= V / dri, ( J;\k,i+T)^p\~k(rt),T) 

L m =l J ° 

xF( im '(fe(f(t + T))). (47) 
4.2 Long time tails 

VACF has already been obtained for the dilute granular 
gases under a simpler approximation in Ref. |13j . but we 
can discuss them by the more precise treatment mentioned 
in this paper. 

To obtain Cs, a a(t) in Eq. (U)), we need to solve the 
time evolution of P(k, t) = P(q, t) given by 

d f -i*kx-^+D*k 2 )p(k,t) = 0, (48) 



with D* = (a 2 tE)D. The solution of this equation is 
given by 

P(k,t) = P{k{rt),0)E D (k,t), 

with 



En{k, t) = exp 



dsD*k{Y{i~ s)f 



(49) 
(50) 



Hence, we obtain the time correlation of P(q, t) as 
[P(q^P(q',0)) = (p(fc,t)P(fc',0)) 



(51) 



where we have used (^P{k, 0)P(fc', 0)y = 5 k i _ k ^*^. Then, 

substituting Su a (q,t) — (t^ 1 a 4: )z a+ i(k : t) into Eq. (0} 
with Eqs. (1ST1) and (|46)l . we obtain 



/dk f dk' 



(2tt) 3 J (2tt) 3 
z a (k, t)z a (k', Q))E D (k, f)£fe',-fe(7*t)' 

/dk 
j^C a+1 a+1 (k, t)E D (k, i),(52) 



Substituting Eq. (gTJ) into Eq. {55]) with the aid of Eqs. 
(HP and (EOl), we obtain 



(53) 



with 



, , drE^(k,t + r)E^(-k(j*^,r) 

Jo 

E D (fc,t). (54) 



It should be noted that Cg ™^(i) for 2, m = 1, 2 is the corre- 
lation of the longitudinal components of the velocity, while 
^s'qq(^) f° r m = 3, 4 is the correlation of the transverse 
components. Cg'^W f° r other cases is zero because p(' m ) 
in Eq. (|54p is zero as shown in Eq. f|45|) . 

It should be noted that Eq. (|53")l describes the veloc- 
ity autocorrelation functions for the sheared granular flu- 
ids for all-time region. In the short-time regime t < j^ 1 
Cs,aa(t) in Eq. (|53")l . is reduced to the known result for the 
fluid at equilibrium Cs, aa (t) ~ t~ z / 2 . For t > 7 -1 , they 
decay faster than t~ 3 / 2 . The detailed calculation of the 
time correlation function in the long-time region is sum- 
marized in Appendix |DJ Here, we only present the results 
of our calculation in this subsection. 

The longitudinal components of the velocity autocor- 
relations Cs™a(t) with l,m — 1,2 approximately satisfy 



= 0{{^' 2 )- 1 ), 



,(22) 



(12) 



.(21) 



(55) 



l.m — 1 



See (fTU7|) - (fTU^j) in Appendix ID.2I for the explicit expres- 
sions of Cg '^(f) with I, m = 1, 2. 

The transverse modes of velocity autocorrelations Cg ™^ (t) 
with I, m — 3, 4 are approximately given by 

(56) 

Oo = 0') = 0') = o (57) 
42(0 « (^ 5/2 )- 1 , 

42(0 « (7* 972 )- 1 , (58) 

The explicit expressions for Eqs. (I55|) ~ ((55)) are given by 
Eqs. ([TT2"|t . dBp in Appendix E3 

Substituting Eqs. (|55 ]) -([55 )l into Eq. ([55|. the long- 
time behaviors of the velocity autocorrelation functions 
are given by 

C s<aa (t) oc (t* 6 / 2 )- 1 . (59) 

We should note that the long-time behavior of Cs, QQ (i) 
is previously obtained in Ref. [13] , where the mixing term 
l_i in Eq. (fTT)) is ignored with the method for the kinetic 
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equations for dilute granular gases relying on |33j . The 
exponent for the long-time tail in Eq. ([59]) is identical to 
that in Ref. [T3], but the amplitude for Cs,aa(t) should 
be different. The validity of the theoretical prediction for 
G's,aa{t) has already been verified by the numerical simu- 
lations with the introduction of one fitting parameter |13j . 
The results of the two-dimensional simulation reasonably 
agree with the theoretical predictions, that is, (i) the more 
precise treatment mentioned in this paper gives us only 
corrections of amplitude of the autocorrelation function, 
and (ii) the previous results can be used for dense gran- 
ular systems. If we ignore the second term on the right 
hand side of (1251) (|27|) which represents the mixing effect, 
Cs,aa(t) is reduced to that in Ref. [13] . 

Our prediction t~ 5 ^ 2 is different from i -9 / 2 [8,9 j. We 
also stress that our theory gives reasonable agreement 
with the simulation including the existence of anistropy in 
VACF and the cross-over regime around t ~ -y -1 for two- 
dimensional case|13). Therefore, we believe that the ve- 
locity autocorrelation function satisfies t~ 5 ^ 2 for t ^> 7 _1 , 
but the direct check of its validity based on simulation for 
three dimensional systems is difficult. 



4.3 long-range correlation function 

Our formulation can be used to discuss the equal-time 
spatial correlation functions. It is obvious that the formu- 
lation of our theory is almost the same as that in Ref. 
[I?] , and the results is identical to those reported in the 
previous paper. Therefore, we only present the result of 
our analysis in this paper. 

Since Sh(q, t) — z\{k, t) and Su a (q, t) = t^, a 4 z a +i[k, t), 
the spatial correlations in Eqs. (|5|) and ([6]) are given by 

dk dk' 



correlations for Eqs. (|6"2"|) and (|63p are given by 

-11/3 



C n „(r) oc 



C PP (r) oc ( - 



-5/3 



r > l c , 
r > lc, 



(64) 
(65) 



which indicate the existence of the algebraic correlation 
in C nn (r) and C pp (r) with l c = a/j*. The validity of the 
existence of long-range correlation has been verified by the 
numerical simulation in Ref. [T7J , where the theory can be 
used, at least, until <fi < 0.50 with the volume fraction 4>. 



C nn {r) = er~ 



Cpp(r) = (mn at E ) 



(2tt) 3 (2tt) 3 

i{k-r-\-(k-\-k')-r' }/ct 

dk dk' 



z 1 (k,0)5z 1 (k',0)) 



(60) 



(2tt) 3 (2tt) 3 
(2«(fc,0)2a(fc',0)) 

a=2,3,4 
— i{k-r / a-\-(k-\-k')-r'} /a 



(61) 



Substituting Eq. gBJ into Eqs. §§§ and dHTJ), we ob- 



tain 



C„„(r) = u~ 



dk 

72^ 



Cn(fc,0)e- 4fcr / CT , (62) 



/dk 
—-C aa (k,Q) 
. V / 



a=2,3,4 * 



-ikr I 'a 



(63) 



where the definition of C a p(k,t) is given by Eq. (|47|) . As 
shown in Ref. [17], using the similar scaling procedure to 
obtain Eq. ([55]) . the asymptotic behaviors of the spatial 



5 Discussion and conclusion 

In this paper, we demonstrate that an unified descrip- 
tion of both equal-time spatial correlation functions and 
the time correlation function for the uniformly sheared 
granular liquids with small inelasticity is possible based 
on the generalized fluctuating hydrodynamics associated 
with FDR. The theory can be used for (i) an uniform 
sheared state with a constant temperature is stable, and 
(ii) FDR and Navier-Stokes like equation can be used. 
These assumptions are justified for nearly elastic granu- 
lar liquids in a small box under Lees-Edwards boundary 
condition. The quantitative validity of our treatment has 
already been confirmed in Refs. [T3] and [T7J. The theo- 
retical prediction for Cs, aa {t) in this paper is essentially 
the same as that in Ref. [13] except for the amplitude. 
Although the method in Ref. [13] is only valid for dilute 
granular gases with small inelasticity, we now obtain more 
generalized method which can be used for considerably 
dense granular gases. Indeed, semi-quantitative validity of 
our method for the equal-time spatial correlation functions 
has been confirmed for tj> < 0.50 [T7J. We also note that 
the method in Ref. [13) ignore the mixing terms among 
transverse velocity mode, the longitudinal mode and the 
density mode, but we now include such mixing effect in 
our calculation. Thus, we conclude that the contribution 
of mixing terms is small. 

Although we have fixed the temperature in our anal- 
ysis, the temperature fluctuates in granular liquids. How- 
ever, if we take into account the fluctuation of the temper- 
ature, the asymptotic behaviors of the correlations does 
not change. Indeed, in Ref. [IB], we have analyzed the 
long-ranged correlation function from the fluctuating hy- 
drodynamics with a local transport law, where the tem- 
perature fluctuation is considered, and demonstrated that 
the asymptotic decay of the correlation is the same as Eq. 
{BSD. 

We have three directions to extend our work. One is 
to remove FDR in Eq. (flU)) to describe highly dissipative 
granular liquids, another is to use soft-spheres to describe 
very dense granular liquids including jamming transition 
and the other is to include the temperature fluctuation as 
well as nonlinear hydrodynamic equations to describe non- 
uniform sheared granular liquids. The last one should be 
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important to extract the characteristics of granular liquids 
under physical boundaries. 

In conclusion, we demonstrate that an unified descrip- 
tion of both time correlation function and spatial corre- 
lation function is possible for uniformly sheared granular 
liquids. This theory unifies our previous work written in 
Refs. [13] and [17]. This method is also valid for sheared 
molecular gases controlled by the Gaussian thermostat. 
This unification may give us a simple view of sheared gran- 
ular liquids. 

We thank S.-H. Chong for discussions. This work was 
supported by the Grant-in- Aid for scientific research from 
the Ministry of Education, Culture, Sports, Science and 
Technology (MEXT) of Japan (Nos. 21015016, 21540384, 
and 21540388), by the Global COE Program "The Next 
Generation of Physics, Spun from Universality and Emer- 
gence" from MEXT of Japan, and in part by the Yukawa 
International Program for Quark-Hadron Sciences at Yukawa 
Institute for Theoretical Physics, Kyoto University. 



where we have used Sp a (q,t) ~ mnoSu a (q,t) with no = 
N/V, (P{-q,i)P(-q',0)) = (P{q, t)P(q', 0)), the decou- 
pling approximation, and replaced q' by — q'. Finally, re- 
placing the sum over q by an integral as J2 q ~ * V J dq/(2n) 
we obtain the approximated expression of Cs aa (t) in Eq. 



B Eigenvectors and eigenvalues 

In this appendix, let us explicitly write eigenvectors and 
eigenvalues for the eigenequation (|24p derived from the 
linearized fluctuating hydrodynamics. 



In order to obtain xp^\k), 
use the expansions 



ip u \k), and \ {j \k), we 









(69) 




= - 




(70) 




= ^ H 




(71) 



A The approximated expression of C s , aa (t) by 
the velocity field 

In this appendix, let us derive the approximated expres- 
sion of Cs,aa(t) given by Eq. The derivation is based 
on the explanation in Ref. [26] . As a result, we can calcu- 
late Cs,aa{t) by an unified method. 

In order to obtain the approximated expression of Cs, a a 
we introduce a complete orthonormal set of functions ipi (t) 
of the phase space point satisfying ((pi(t)^(t)) — 5i, n . Us- 
ing 4>i(t), Cs, aa (t) in Eq. ([T]) can be represented as 

C S , aa (t) = ^2 (<5«l,a (0) <Pt (0) ) (</> m (0)6vi <a {0))((pl (t)<Pm ) , 



in terms of 7*. 

We should note that the perturbation in terms of 7* 
is not the expansion from an unsheared state of granular 
liquids. Indeed, it is well-known that properties of sheared 
granular liquids completely differ from those of freely cool- 
ing granular liquids. In the case of sheared granular liq- 

2 from 



uids, we obtain the relation 7* ~ 7/VT 1 ~ — e 
the balance between the viscous heating and the collisional 
.energy loss. Thus, the expansion in terms of 7* can be re- 
* garded as that by small inelasticity [32] . 

Substituting Eqs. ([69]), ([70]) and ([71]) into Eq. ([24]), we 
obtain the zeroth and the first order perturbations as 



(Lo 



A«l).^(fc) 



0, 



(L -A£>l)-^>(fc) 



(66) 

where we have used Yli(^ v i,ct^ v i,a) = N(6vi ja 5vi tCt ). 

To extract the slow relaxation of Cs, Qa (t), we choose 
the slow variables 



-lk x 



dk„ 



Li-A^l 



(72) 



^\k)=0. (73) 



.(*) 



VmJVT 



5p a (q,t)P(-q,t), 



(67) 



Solving these equations, we obtain the eigenvalues given 
by Eqs. ([25]) - ([28]) . Similarly, we obtain the right eigenvec- 
tors 



(1)T 



■ / . , 3 , k x ky k z 
— — iknocr , A+ — , A+-p, A+ — 



where 6p a (q,t) = m J2 3 ■ Svj t0l (t)e "J-*"^*) and P{q,t) — 
e -ig-ri(i) are |- ne p 0UI .j er components of the momentum 
density and the microscopic concentration of the tagged 
particle, respectively. We should note that (5v\^ a (j)* q p) = 

y/T/(mN)6 a ,p. Then, substituting Eq. (JBTJ) into Eq. (jgrJJ). 
we obtain 

C s , aa (t) ~ ^(^ 1 , Q (O)C(g,O))(0 Q ( 9 , ,O)^ 1 , Q (O)) 
(0 Q (q,t)C(g',O)), 

- y^^2( 6{i 'o l (q,t)Su a (q' ,0)) 
q,<i' 

(P(q,t)P(q',0)}, (68) 



1> 



(2)T 



1 

1 

iv_ 

#(3) 
3/(4) 



— — iknoa ,A_ — ,A_-p,A_ — 



M{k)& 



(4) 



(1) _ 



and the left eigenvectors 
1 

1 

nZ 

tp<& = <p( 3 ), 

<^ 4 > = -M(fc)# (3) + 3> (4) 



ikp*(k, e), A+ 



k x 

k x 

T 



,+ k ' 

A_^,A„ 



k 

kz 
' k 



(74) 

(75) 

(76) 
(77) 

(78) 

(79) 

(80) 
(81) 
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where we have introduced 

gf(3)T _ <p(3) _ f g k± k y k z 

kk± ' k ' kk± 

k i ki 



and 



M(fc) 



■tan 1 (k y /k±) 



(82) 
(83) 

(84) 



with k± = k 2 — ky. It should be noted that these eigen- 
vectors with M(k) is valid only when k x ^ 25 . 

C Solution of Eq. flU 

In this appendix, we derive the solution (|38l) of Eq. (l36|) . 
Introducing Q = k(j*F) and transformations (t, fe) = 
(F, Q{-i*t')) with P = f and Q(f) = (Q B , Q y + f Q x , Q 2 ), 
Eq. (1361) is rewritten as 

(c%, + AW(Q(-7*?)))) A^{QJ) = F^(Q(~rt')),i'), 

(85) 

where we have introduced A^(Q,F) = a& (Q(--y*F), F), 
and used the relations c% = c%/ + j*Q x dQ y . 
The solution of Eq. ([55} is obtained as 

E'U\Q,F) 



A^(Q,F) = / dflFWCOC-r*),*)^ 

A' E 



+ ^(Q^)^ )(Q ' F) 



'^(Q,?o)' (86) 

where is an initial time, and E ^'(Q,t) is defined by 

E'^(Q,F) = exp[- [ dfA«(Q(-7*f))], (87) 
Jo 

Taking the limit F Q — > -co, we obtain 

i«...(0. ? , = /'j^.(0(- 7 .»,,»,|^|fl. 

Since a«(fc,t) = A^(k(j*t),F), 



£'«(fe(7*t),s) 
Here, from the transformation r' = t — r, we obtain 
E'W(fe(7*f),t) _ exp[-/ *dfAW)(fc(7»(t-r)))] 



E'^{k(j*F), s) exp[- J S dfA0')(fe( 7 *(t - f )))] 

= exp[-/Jrff'A»>(fc( 7 V))] 
exp[-/^ s df'A0")(fe( 7 *f'))] 

pi— s 

= exp[- / df'A^')(fc(7*f'))] 



(90) 



where E^(k,i) is given by Eq. ([391) . Substituting this 
equation into Eq. (j89[) . we obtain Eq. (j38[) . 



D Long-time behaviors of the time 
correlation functions 



In this appendix, we explicitly evaluate the long-time be- 
haviors of the time correlation functions given in Sec. 14.21 
In the first part, we obtain the general expressions for the 
long-time behaviors of the time correlation functions. In 
the second part, we evaluate the long-time behaviors of 
the time correlations of the longitudinal component of the 
velocity. In the third part, the time correlations of the 
transverse component of the velocity is estimated. 



D.l The general expressions for the long-time 
behaviors of the time correlation functions 



Introducing K as K x = k x j*P^ 2 and K a = k a FI 2 for 
a/i, and t' = r/t, Eq. ([54l) are replaced by 



dr' 



W J 

F {lm) {K + F- 1/2 (1 + t')K x ) 

(K, t(l + r'))£ (m) (-K - r 1/2 e y K x ,tT') 
^ +l (k)^[(-k-t-" 2 e y K x ) 
E D (K,t). (91) 

where we have introduced K = t~ 1 / 2 (K x /(j*t),K y ,K z ). 

In the long-time regime t ^> 7 _1 , since k — > 0, E"'(k, ir') 
and Ed(K ,tr') are approximately given by 



\E^\k,i r ') 



exp 



ds'vl{0,e)K T {s'f 



K t (t 



\E^(k,ir') = \E^(k,tr% 



7T'( 92 ) 
(93) 



E^{k,ir') 



exp 



'\2 



ds'u^(0,e)K T (s') 



E (i \k,ir')^E^(k,ir') Kt 



K t (t') 
K t 



E D {K,F)~exp 



K T {T>y 

ds'D*K T (s') 2 



(94) 



(95) 



where we have introduced K t (t) — (K y + tK x ) 2 + K 2 
and Kj, implies Kt(0)- 
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ip"'(K) in the long-time regime is approximated by where we have introduced 



V>« T (K)~-±={#o+^}, 
^ 2 > T (K)~-L{#o-^}, 

</> (3)t (k) 



(96) 
(97) 



0, -M'(K) 



K z \K Z \ K y K z 



\K z y K T ' K T \K Z 



^ T (x)~(o,- T ^ I1 o ) . K ' 



\Kz\ i*t\K z 



with 



^o- (1,0,0,0), 



K x K y K z 



(99) 

(100) 
(101) 



where we have introduced M'(.K') = ^jp^-j tan 1 
Here, we have introduced the infrared cut off A for fc^ to 
avoid the divergence of M'{K). 

p(im) ^ j n ^j^g long-time regime is approximated by 



^(22) = £,(12) = ^(21) = 



(102) 

(103) 



(104) 

F^iK)- -N(K)F^ 3 \K), (105) 
F<M(K) = -F&Xk) ~ M'{K)F^\K) (106) 

with iV(K) = M'{K) 2 + 1. 

D.2 The time correlations for the longitudinal 
components of the velocity fields 

Substituting Eqs. p ]l - lfT06)l into Eq. fEE]), <?s™2(*) for 
/, m = 1, 2 in the long-time regime is approximately given 
by 

ICffiWI < ^f^(7-? /2 )- 

J l o °° dr'G^(K,r'), (107) 

and 

i<Lwi < ^M^r 1 

j-^j~dr'GW(K,T') (108) 

for q = y, z, and 

i4"L(*)i = ic£L(*)i = = iqgLwi, 

(109) 



G ^ (JC, = H X (K, r>) K f + ^ , (110) 



with 



Hi {K, t') = exp 



2 Jo 



l+r' 



'\2 



exp 



exp 



dsV*(0,e)K T (s') 



ds'vl{Q,e)K T {s'f 



[ ds'D*K T (s" 2 
Jo 



(111) 



Eqs. (|107l) - (|109p are the explicit expressions for the long- 
time behaviors given by Eq. (|55|) in Sec. 14.21 



D.3 The time correlations for the transverse 
components of the velocity fields 

Substituting Eqs. @-(HM]) into Eq. flU]), G^(t) for 
I, m = 3, 4 with a — x, y is given by 

li^C dT ' G « l){K ^ (112) 

where we have introduced 
G' 33 ) (Jf, 7-0 = ff 2 (A-, r0 ^(l^+r0 4 T( ^ Q)T(K) 1); 

(113) 

Gf)(K,rO =g2 (K,rO ^; ( y- , (114) 



r') 2 + 1 



H 2 (K,t')K t (1+t') 2 , (115) 
G^\K,t>) ~ - ^ r(1 A + T/)4 T(A:,l + rOT(A:,0) 



(116) 



Gi 43 > (K, /) ^ _ ^(^+ r/ ) r (JC| x + / )T(K) 1} 

H 2 (K,t'), (117) 
Gf)(K,rO = G( 34 )(K,rO = Gf^rO = 0, (118) 



with 



H 2 (K,r') = exp 



dsV*(0,e)if T (s0 2 



exp 



exp 



>V2 



ds'v%(0,e)K T (s') 



[ ds'D*K T {s'f 
Jo 



(119) 
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and 



(120) 



Eq. (|112l) is the explicit expression for long-time behaviors 
given by Eqs. ([56]) and dSZ} in Sec. PI 



Similarly, from Eqs. ^-fMfr into Eq. G^jO) 
for Z,m = 3,4 in the long-time regime is approximately 
given by 

c m (t) = A^ ( ., 3?/2rl 

C (. 2(t) = i^ ( ., 2F/2r , 

/ (i2i » 

where we have introduced 

Gi '(JSf.rJ-flbCJC.r) ^ r(1)2 . 



(122) 

Gi 33 ) (Jf, r') = | ^±^1!t(K, 1 + r') 2 + 1 



(JT, r') = H 2 (K, r') KTi \ + ^ vKx T(K, 1 + r'), 

(124) 

(125) 

Equation (|121l) is the explicit expression for long-time be- 
haviors given by Eq. (|S"8|> in Sec. 
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